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Abstract 

Layered and nanotubular systems that are metahic or graphitic are known to exhibit unusual 
dispersive van der Waals (vdW) power laws under some circumstances. In this paper we inves- 
tigate the vdW power laws of bulk and finite layered systems and their interactions with other 
layered systems and atoms in the electromagnetically non-retarded case. The investigation reveals 
substantial difference between 'cleavage' and 'exfoliation' of graphite and metals where cleavage 
obeys a C2D^^ vdW power law while exfoliation obeys a C3 log(L'/Do)-C'^^ law for graphitics and 
a 6*5/2-0^^/^ law for layered metals. This leads to questions of relevance in the interpretation of 
experimental results for these systems which have previously assumed more trival differences. Fur- 
thermore we gather further insight into the effect of scale on the vdW power laws of systems that 
simultaneously exhibit macroscopic and nanoscopic dimensions. We show that, for metallic and 
graphitic layered systems, the known "unusual" power laws can be reduced to standard or near 
standard power laws when the effective scale of one or more dimension is changed. This allows 
better identification of the systems for which the commonly employed 'sum of CqD~^^ type vdW 
methods might be valid such as layered bulk to layered bulk and layered bulk to atom. 
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Layered bulk systems such as graphite and boron nitride have their atoms confined to 
a series of spatially discrete planes with interplanar distances significantly greater than the 
intraplanar atomic separations eg. 3.34A vs 1.41A for graphite. They can exhibit unusual 
electronic behaviour due to the nanometer scale of the layer thickness and macroscopic 
scale of the other two dimensions. This scale variation means that great care must be 
taken in investigating subtle physical effects such as dispersion forces as previous examples 
demonstrate 

All separated electronic systems exhibit long-range attractive potentials arising from in- 
stantenous electron fluctuations correlating via the Coloumb potential. These long-range 
potentials (often called van der Waals dispersion potentials when electromagnetic retarda- 
tion is ignored) are typically absent from the commonest ab initio calculations such as DFT 
in the LDA or GGA, or are approximated by pair-wise inter-atomic potentials of the form 
CqD^^ which are 'summed over' in some way (see eg. refs 0, [s], [g], [g, [?, ^, ol) to obtain a 
new effective power law of the form C„i5^". Here the exponent n is an integer and depends 
only on the geometry of the system while C„ depends on the individual atoms as well as the 
geometry. ^ 

As summarised in [1] the exponent of an asymptotic power law in metallic and graphitic 
systems can depend on both the geometry and the type of material, with metals, graphene 
and insulators all differing. For example, with two parallel, nano-thin layers of a metal, 
graphene and insulator, the power-law exponents are 5/2, 3 and 4 respectively (where insu- 
lators do obey a 'sum over D~^^ rule). When the number of layers is inflnite we will show 
that the power depends not only on the type of material, but also the way the layers are 
divided. This will be investigated through three types of division: equal separation of all 
layers ('stretching'), division into two sub-bulks ('cleavage') and removal of one layer from 
the top of a bulk ('exfoliation'). Furthermore the interaction of layered bulks with atoms 
will be studied. 

In this paper we investigate metallic and graphitic systems under these different types of 
division (insulators have trivial 'sum over atoms' exponents and need no further investiga- 
tion). Neither metals nor graphitic systems are guaranteed to obey 'sum over layer' power 
laws and special care must be taken to evaluate their long-range correlation effects. As with 

previous work , Q, Q, Q, Q, Q, H, Q, 0, Q we make use of the Adiabatic-Connection 

Formula and Fluctuation Dissipation Theorem (ACFFDT) under the random-phase approx- 
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imation (RPA) to calculate the leadmer power laws under these different methods of division. 

n 

All results in this paper are for the electromagnetically non-retarded case which [19| show 
to be unimportant in the range of interest for similar systems. 

Stretched graphitic systems ('straphite') have already been studied in jsl where it was 
shown that the dispersion potential for an infinite number of graphene layers, each with inter- 
layer spacing D follows a CsD~^ asymptotic power law at T = OK where C3 = O.SOeVA . We 
may make use of the same basic approach employed for graphitic layered systems to calculate 
the dispersion for metallic layered systems (we believe that graphite-metal intercalates may 
be examples of this type of system). For brevity we define an intralayer Coloumb potential 
multiplied by a density-density response function 

C{q,u) =x{q,u)w{q) (1) 

where q is the wavenumber parallel to the plane, xiq-,'^) is the non-interacting electron 
density-density response in a layer and w{q) is the Coloumb potential. When the system is 
metalhc this takes the form C{q,u) = ~uj^2D{qY/u^ for g — > where uj^2T>{q) = { 2egrn )^'^^' 
riQ = Ai'e- layer Miayer is the 2D electrou density of each layer and rUe is the mass of an 
electron. We may thus utilise equations 4 and 7 of ref. and change variables to = qD 
and T] = \/Dujp2Diq)^^u to show that the difference in the correlation energy per layer from 
the infinitely separated case takes the form 

1 /»1 roo roc 

f^vdw-^p^y^ dAy^ O'^'dej^ dr]gie,r],X) (2) 

demonstrating a C5/2-D"^^^ dispersion power law for layered metals. Numerical evaluation 
gives C5/2 = 9.26/l( ^"°^ Y^^. As with graphitic systems this power law (although not the 
constant prefactor) is universal to all multi-layered metallic systems with an infinite stack 
of isotropically stretched layers. 

Given the universality of the van der Waals exponent for isotropic stretching of a given 
material it is worth exploring the validity of a 'sum over rule for layered systems. In 

Table [I we present the ratio of the C„ coefficient for a tri-layered or stretched system to the 
bi-layered system for metals and straphite as well as a the 'sum over C„' prediction for this 
ratio (here n is 3 for straphite and 5/2 for metals). The prediction is somewhat sound for 
straphite with an overprediction of 14% for straphite but is much less so for metals where 
it leads to a 30% overprediction of the potential. This suggests that, even using a correct 
power law, rules which effectively sum the coefficients may prove troublesome. 
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System Cn'^'^*''™/C^'Predicted Error 



Tri-graphene 


1.3632 


1.4167 


3.6% 


Stretched 


2.1157 


2.4041 


14% 


graphite 








Tri-metal 


1.3496 


1.4512 


7.5% 


Stretched metal 


2.0628 


2.6830 


30% 



TABLE I: Comparison of the system to bi-layer ratio of the per-layer vdW coefficients for graphitic 
and metallic systems, for tri-layered and stretched systems. Column 3 is obtained by summing 
over the individual layers and takes the form |(2 + 2^"') for tri/bi and 2({p) for stretched/bi. 




FIG. 1: Layers after 'cleavage' 

'Cleavage' represents another means of division of a layered system. Here the system is 
split between a single pair of layers to form two new layered systems. We refer to the new 
systems as half-bulks as opposed to the original full bulk. For homogenous, infinite layered 
systems the two are mirror images of each other. 

Separating a layered bulk into two smaller half-bulks is equivalent to keeping all but one 
layer at an inter-layer spacing d while increasing the remaining one to D ^ d as in Figure [H 
(we ignore any relaxation of the layers or layer spacing at the newly created surfaces). Here 
we are interested in the total dispersion energy per unit area rather than that per layer per 
area as in the previous case. 

We may make use of the second-order perturbation formula for the dispersion energy 
between the half-bulks while treating that within them to all orders (sometimes called the 
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Zaremba-Kohn 20] formula). Thus 



oo 



f^vdw = v4 / / qdqT{q, u; D) (3) 



where 

^(g, n; Z^) ^ e-^'^^Tr [{Xx=i{q, rf)^'(g^))'] (4) 

where w\j = we~'^^^^~^^^ governs the interaction only between the two separate half-bulks 
while Xij is the full interacting response of the density in layer i to a potential perturbation 
at layer j within a single half-bulk. 
Expanding the trace gives 

ijkl i 

where Sui is the interacting electron density response of layer i in the half-bulk to an external 
potential perturbation 5vi = wp^ with p = e~'^'^. 

We may calculate Sui through the RPA equation Sui = Cp^ + CJ2j>oP^^~''^^''T'j where 
C = C{q,u) is as defined in equation ([T]). This gives the following recursion relationship 

+ Sui = [{p + p"^) - C{p - p"^)] Sui+i (6) 

where 6nQ = C(l + J2j>oP'^^^j) require Suoo < oo. Writing a formal power series 

generating function Jli{x;C) = J2i>o^''^i^^ transforms equation into 

(1 + x'^)Ai{x] C) — — x6ni = 

X [{p + P~^)+C{p-p-^)] {M{x-X)-6no) (7) 

so that (with correct series asymptotics) 

Mix;C)= ^^^'f'^ (8) 

p[l — r_x)[l — r^p ^) 



where r± = ( ± a/C^ — 1, C = \{p + P + \^{P ~ P Noting that ^^p^Sui = M.{p) lets 
us write equation ([5]) as 

T{q, u; d, D) = e-^'^^'Mie-'^^; C{q, u)f. (9) 
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For graphite shows that C{q,u) = k[1 + u'^/{voq)'^]~^^'^ where Vo = 5.0 x lO^ms"^, 
K = j^j— = 12.1. Defining a = d/ D, 9 = qD and r] = ud/vo lets us expand T{9,7]) = 
J-'{qD, ud/vo] a, 1) in powers of a. To leading order this gives 

HO, V) = e-''- "l^^,, (10) 

for graphite. A similar change of variables yields a similar leading order expansion for metals. 
For these cases equation ([3]) gives a leading power law of exponent = 2 as with insulators. 
Thus for grapite and metals we can write 

U,^^{D ^ oo) ^ ^ (11) 

where 

c^^ = iJ^Pf d,r(„). (12) 

HereT(r7) = ^k^j^k + r/ + a///^ + 2Kr/)-2 for graphite and T(r7) = ( 1 + + a///^ + 2)^^ 
for metals with cUpSD = (^^^f^)^*^^- graphite we find = 0.13d~^eVA^ and for metals 
C2 = 1.6886 X 10^^fkjp3B- The latter result agrees exactly with continuous but anisotropic 
models of half-bulk metals where electron movement is restricted to be parallel to the surface, 
as is expected if the d ^ limit is well defined. 

These results differ significantly from those expected by a simple sum-of-layers approach 
where we would expect graphene to obey a CiD~^ power law, and metals to obey a Ci/2-D^^/^ 
power law. The screening in these layered systems seems to cancel the different correlation 
effects of the individual layers so that they act as pseudo-insulating bulks. 

Another means of dividing graphene ('exfoliation') is to peel a single layer of graphene 
from the top of a bulk. This represents yet another method of division where one system 
is a layered half-bulk and the other a single layer. We restrict our investigation to the case 
where the removed layer plane is always stiff and parallel to the planar surface of the bulk. 

To model exfoliation we use a similar perturbative approach to cleavage but in equation 
([3]) set = e~'^'^^Ai{e^'^^;C)j^ the product of the interacting response of a half-bulk and 
a single layer. For graphene this complicates the problem as the a = d/ D ^ limit works 
well for Ai but not for Here we keep a in the single-layer response only. Under 

transformation of variables 6 = qD and 1] = (which already gives a outside the 
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integrals) we find 



X (V?7^ + {aey + Kaey\ (13) 



Setting 9 = 1 (where 6'^e~^^ takes its maximum), approximating ^Jif + by a + 77 and 
K + T] + \/ 'ff + 2^1] by 2{k + 77) allows us to approximate the rj integral to show a leading 
a log(a;) term. This is equivalent to a power law of the form 

U^dw[D ^ 00) ^ — . (14) 

and numerical calculation of equation ([3]) validates this assumption. For graphene (where 
K = 12.1) we find Do = O.lGd and C3 = O.OTeVAl 

A similar analysis of metals shows a C5/2-D~^/^ power law with C5/2 = 6.42 x 
lO^^^/dhwpSD = 9.07 X 10^^hujp2T){(l = 1). The 5/2 power law is the same as that of a 
metal layer interacting with a continuum model of a metallic half-bulk as is, again, pre- 
dicted by the limit of 0. 

While this analysis is not valid in the small D regime it is worth noting that exfoliation 
and cleavage exhibit different power laws for graphite. This suggests that sum of models 
or converting experimental results from one to the other such as those employed in ref. 



2l| may need reexamination. Unfortunately accurate calculation of the dispersion energy 
of such systems for D ^ do (where do is the layer spacing of graphite) is as yet intractable. 

In nanoscale systems there are often combinations of molecules, layers and bulks. A 
simple example is an atom interacting with the surface of a layered metal or a molecule 
interacting with a graphene surface. Here the power law could be affected by the layering 
and electronic properties of the material. 

In the coordinate system used for the layered models the interacting response function 
of an infinitely small "atom" located at Zz can be written as 

XA(q,q',2,^';M) = 

a{iu) [q ■ q' + d,d,>] 6{z - Z)6{z' - Z) (15) 

where q and q' are reciprocal lattice vectors in the plane and a{iu) is the interacting dipole 
polarisibility of the atom at imaginary frequency u. This formula, used in equations 
correctly reproduces the CqD~^ power law for interacting atoms. 



Graphite Metal 



Insulator 



Stretching D''^ D'^/'^ 
Cleavage D'"^ 
Exfoliation log{^)D^''^ £)-5/2 
Atom-bulk D'^ D'^ 



^-4 



TABLE II: Asymptotic power laws for various systems demonstrating both the material and struc- 
tural dependence. 

We can use equation f|T5l) to calculate the interaction of an atom with a layered bulk 
(metallic or graphitic) by making use of equations ([3lll]). This gives a power-law exponent 
n = 3 in agreement with the prediction of a sum over CqD~^ potentials. This result strongly 
suggests that the unusual power laws exhibited by layered systems result from the interaction 
between long-range fluctuations in both systems and that removing them from one reduces 
the systems to 'typical' dispersive behaviour. 

As has been seen here and in other work the asymptotic power law behaviour of 

layered systems can be anything but simple. Both graphitic and metallic systems exhibit 
vastly different dispersive power laws to insulators so that 'sum of Cg' approximations such 
as those typically employed cannot be used in the asymptotic region. It seems unlikely that, 
with such varied asymptotes, the cohesive energies and other similar measurables can be 
investigated using simple models. 

To illustrate these discrepancies we present in Table [III a summary of the various power 
laws studied here. The insulator result represents the 'classic' sum over atomic power- 
laws behaviour of each system and any difference from its exponent represents 'unusual' 
behaviour. 

While uniform stretching results in different power laws for metals, graphitics and in- 
sulators, cleavage removes this variation and involves the same exponent for all materials. 
This suggests that the interlayer screening induced by the Coulomb potential dominates the 
local response of each layer in the van der Waals energy for such systems converting their 
behaviour into that of non-layered or insulating bulks. This is further demonstrated by the 
fact that both the power law and coefficient of cleaved layered metals is the same as that of 
cleaved bulk metals with electron movement restricted to the plane. 
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By contrast, keeping a finite number of layers asymptotically isolated, as in exfoliation, or 
all layers asymptotically separated, as in stretching, returns different power-laws for different 
systems. In these cases at least one layer can be considered infinitesimally thin which we 
believe to be a requirement for the unusual power-laws. Replacing the isolated layer by an 
atom, however, returns the classical results which suggests that at least one large dimension 
is required for unusual vdW dispersionpower-laws as postulated in refs. l2| 

Overall, as this work and references demonstrate, the dispersion forces of systems 

with a mix of nanometre and macroscopic length scales are more complex than classic Lifshitz 
theory predicts. As the differing power laws for cleavage and exfoliation demonstrate we 
must take great care in using indirectly derived cohesive energies from experiment. 

These unusual van der Waals power laws may also have profound effects on the behaviour 
of many nanosystems. For certain systems it may be neccessary to adapt molecular dynamics 
and other semi-empirical and approximate ab initio simulation methods to account for these 
differences in order to best replicate experiment. 
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